oo 



a Department of Physics, 



BU-HEPP-03/12, UTHEP-03-1201 

Dec. 2003 

Comparisons of Fully Differential Exact Results for O(a) Virtual 
Corrections to Single Hard Bremsstrahlung in e + e~ Annihilation 

at High Energies^ 

O ■ C.Glosser abc , S. Jadach de , B.F.L. Ward a and S.A. Yost a 

(N 

3 ! Baylor University, Waco, Texas 76798-7316, USA 

b Department of Physics and Astronomy, 
The University of Tennessee, Knoxville, Tennessee 37996-1200, USA 

c Department of Physics, 
Southern Illinois University, Edwardsville, IL 62026-1654, USA 
OO \ d CERN, Theory Division, CH-1211 Geneva 23, Switzerland, 

0\ ■ e Institute of Nuclear Physics, ul. Radzikowskiego 152, Krakow, Poland, 

O 
^1- 

-»^ ■ Abstract 

We present comparisons of the fully differential exact virtual correction to the important 

qj' single hard bremsstrahlung process in e + e~ annihilation at high energies , which is essen- 

43 tial for precision studies of the Standard Model from 1 GeV to 1 TeV, as calculated by 

> . two completely independent methods and groups. We show that the two sets of results 

are in excellent agreement. Phenomeno logical implications are discussed. 
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Now that the Standard Model electroweak theory has been established at the one-loop 
level [1,2], the stage is set for studying its consequences as both signal and background 
to the physics objectives of current and planned high energy colliding beam devices, from 
1 GeV to 1 TeV for the cms energy in e + e~ annihilations for example. The attendant 
per mille level studies necessitate control of the EW higher order radiative corrections at 
least to order 0(a 3 L 3 ) for the leading log effects and to the exact 0(a 2 ). One set of the 
important contributions to the latter exact results are the virtual corrections to the single 
hard bremsstrahlung in e + e~ annihilations, which therefore have been studied by several 
groups [3-6]. Comparisons of these results are essential in order to gain some confidence 
in their correctness and some facility in their use to confront the SM with precision data. 
For example, in the effort to exploit the radiative return from cms energies in the 1-2 GeV 
regime to the resonance regime of the im system in the Daphne environment, precision 
predictions of the type compared in this paper are essential. Similarly, in order to use the 
radiative return from 200 GeV to the Z in the final LEPII data analysis for precision EW 
tests, again precision predictions of the type compared in the following for the respective 
return processes are essential. 

Indeed, in Refs. [5], some of us (S.J., B.F.L.W., S.A.Y.) have presented comparisons 
of the results in Refs. [3-5] and in general a very good agreement was found. However, 
if one looks at the comparisons in Ref. [5], one can see that at the level of the NNLL 
(next-to-next leading log), there was a difference in the results that was consistent with 
the different levels of "exactness" in the calculations. Specifically, the mass corrections 
are included in Ref. [5] in a fully differential way, whereas in Ref. [3] the mass corrections 
are included but the photon angular variable is integrated over and in Ref. [4] the results 
are fully differential but the mass corrections are incomplete. These comparisons therefore 
can not really test the NNLL, fully differential results in Ref. [5]. 

The situation has changed recently with the advent for the results of Ref. [6]. Us- 
ing a completely independent method and calculation, the authors in Ref. [6] have also 
achieved a fully differential, exact 0(a 2 ) result for the virtual correction to single hard 
bremsstrahlung in the initial state radiation in high energy e + e~ annihilation, with par- 
ticular emphasis on the 1-2 GeV cms energy regime. This then affords a detail cross check 
at the NNLL level of the corresponding results. These comparisons are the subject of this 
paper. 

We note that, already in Ref. [7], a preliminary indirect comparison of the results in 
Refs. [5,6] has been reported via comparison of the two Monte Carlo's PHOKHARA [6] 
and /C/C MC [8] , as these two Monte Carlo's have the realizations of the results in Ref. [6] 
(PHOKHARA) and Ref. [5] ( JCIC MC). Agreement at the per mille level was found on 
selected observables. This is a good basis upon which to view the results which follow. 

Specifically, the Feynman graphs under discussion are illustrated in Fig. ^ I n Ref- [5] 
the respective ISR matrix element is evaluated using the CALKUL/Xu et a/./Kleiss- 
Stirling [9-11] method for the attendant helicity amplitudes using FORM [12] techniques. 
The mass corrections are then added following the methods in Ref. [13], after checking 
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Figure 1: Feynman graphs for the virtual 0(a) correction to the process e + e" 
are illustrated, where f ^ e. 
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that the exact expression for the mass corrections differs from the result obtained by the 
latter methods by terms which vanish as m 2 e / s — > where m e is the electron mass and s 
is the squared cms energy. In Refs. [6], the same ISR matrix element's interference with 
the Born process is calculated in terms of Lorentz covariants in such a way that the mass 
effects are treated exactly. Thus, comparison of the two sets of results gives important 
information on the two different methods of calculation and on the two different treatments 
of the mass corrections. 

For our studies, we will follow the development in Ref. [5] and systematically isolate, 
using the differential cross section , first, the the complete cross section, second, the size 
of the LL contributions, third, the size of the NLL contributions and finally the size of 
the NNLL contributions. In this way, we shall see how well the two independent fully 
differential calculations with completely different approaches agree. This of course will 
reinforce the confidence, whatever it may be, that we have in both calculations. 

Specifically, for the process in Fig. ^ where pifa) is the four-momentum of the 
incoming e + (e~), respectively, and k is that of the emitted real photon, we focus on the 
corresponding virtual correction to the single real bremsstrahlung differential cross section 
da l jdz where da 1 is defined in eq.(3.7) of Ref. [5] and where z = s'/s = 1 — v 
when s is the squared cms total energy and s' is the squared final state // system rest 
mass. Here, we will average over the initial fermion spins and sum over the final ones. 
This cross section has been computed in both Refs. [5] and in Refs. [6] when we restrict 
ourselves, as we do here, to the fully ISR (initial state radiation) corrections and processes. 
Continuing to follow the work in Refs. [5], we again note the that IR (infrared) limit of 
the process is known and can be determined from the arguments of Yennie, Frautschi and 
Suura in Ref. [14]. Thus, with an eye toward the use of the results in Refs. [5] in the YFS 
based CEEX (coherent exclusive exponentiated) Monte Carlo KK MC [8], we work with 
the cross section associated with the virtual correction to the hard photon residual J3[ 
as it is defined in Refs. [8,14] for example. This allows us to focus on the non-IR singular 
part of the respective cross sections. 

What we find is illustrated in Figs. 2-5, for the case // = yT ji + . In Fig. |2J we show 

— (2) 

the complete f3\ distribution for our exact result and our NLL approximate results as 



presented in Ref. [5], the result of Igarashi and Nakazawa et al. [4], the result of Berends et 
al. [3], and, what is new here, the exact result of Kuhn and Rodrigo in Ref. [6]. What we 
see is that there is a very good general agreement between all of these results. To better 
assess the difference between them, we plot in Fig. |3]the difference between the respective 
0(a 2 ) and O{o 1 ) results. Again we see very good agreement between the results. 

To isolate the respective predictions for the NLL effect, we plot in Fig.HJthe respective 
differences between our LL 0(a 2 ) result from Ref. [5] and the other five results. The 
results are still essentially indistinguishable at this level, except that NLL effects become 
apparent in the last data point, for v = 0.9875. 

These comparisons are more evident in the following graphs, Fig. |3] and Fig. EJ where 
we isolate the size of the four NNLL results by subtracting our NLL result from Ref. [5] 
from each of the results. Fig. |H1 is identical to Fig. 03 except for the scale, which permits 
a closer comparison of the NLL results below a cut of 0.95, while omitting some off-scale 
data points beyond that cut. It was already established in Ref. [5] that the expressions 
of Ref. [4] and Ref. [3] agree analytically with our NLL upon taking a collinear photon 
limit. This has been checked as well for the massless limit of the result of Ref. [6]. 

In the NNLL comparisons, all of the results agree to within 0.4 x 10 -5 for cuts below 
0.75. For cuts between 0.75 and .95, the results agree to within 0.5 x 10~ 5 , if the result 
of Ref. [3] is not included, and within 1.1 x 10~ 5 if that result is included. For the last 
data point, at v — 0.9875, the result of Ref. [3] is approximately 1 x 10 -4 greater than 
the others, while the remaining results agree to 3 x 10" 5 . 

These results are consistent with a total precision tag of 1.5 x 10~ 5 for our 0(a 2 ) 

— (2) 

correction [3\ for an energy cut below v = 0.95. The NLL effect alone is adequate to 
within 1.5 x 10~ 5 for cuts below 0.95. The NLL effect has already been implemented in 
the IOC MC in Ref. [8] and the attendant version of KK MC will be available in the near 
future [15]. 

We have also made the analogous study to Figs. 2-5 for 500 GeV. We find very similar 
results, with the total precision tag of 2 x 10~ 4 . 

In this paper, we have presented new comparisons of exact results for the virtual 
correction to the process e + e~ — > // + 7 for the ISR. The results from Ref. [5] are already 
in use in the /C/C MC in Ref. [8] in connection with the final LEP2 data analysis. 

We have compared our results with those in Refs. [3,4,6] and in general we find very 
good agreement, both at 200 GeV and at 500 GeV. The comparisons with the results 
in Ref. [6], which like our results in Ref. [5] are exact but which have the complete 
mass corrections included explicitly whereas our mass corrections are included using the 
approach of Ref. [13] ( which we have shown to differ from the exact mass results by 
terms which vanish as m 2 e / s — > ), allow us to lower our precision tag to 1.5 x 10~ 5 for an 
energy cut below 0.95 compared to what we quoted in Ref. [5]. For example, the size of 
the NNLL correction is now shown to be at or below the level of 1.5 x 10" 5 for all values 



of the energy cut parameter up to 0.95. 

We need to stress that a considerable amount of simplification of large cancelling 
terms in the expressions in Refs. [5, 6] was required to produce the results in this paper. 
Only after all such cancellations were found and carried out analytically could the agree- 
ment shown in Figs. 2-6 be realized. The details of this part of the analysis will appear 
elsewhere [15]. 

Our results are fully differential and are therefore ideally suited for MC event generator 
implementation. This has been done in the /C/C MC in Ref. [8]. The results from Ref. [6] 
are also fully differential with mass corrections and have also been fully implemented 
into a Monte Carlo event generator [6], albeit one without YFS exponentiation. It is 
therefore natural to compare the two respective event generators in the context of the 
radiative return at <J> and B-factories. As we noted above, a preliminary version of such 
results already appeared in Ref. [7], where per mille level agreement was found. The 
results above indicate that, in principle, much better agreement is possible. Further such 
comparisons will appear elsewhere [15]. 

What one sees from the comparisons above is that we now have a firm handle on the 
precision tag for an important part of the complete 0(a 2 ) corrections to the 2/ production 
process needed for precision studies of such processes in the final LEP2 data analysis, in 
the radiative return at $ and B-Factories, and in the future TESLA/LC physics. 
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Figure 2: This is the f3[ 2) distribution for the YFS3ff MC (YFS3ff is the EEX3 matrix element 
option of the /C/C MC in Ref. [8]), as a function of energy cut v m3iX . It is divided by the Born 
cross-section. The exact and NLL results from Ref. [5] are compared to the IN result from 
Ref. [4], the BVNB result from Ref. [3], and the KR result from Ref. [6]. 
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Figure 3: Difference p[ 2) - (3^ for the YFS3ff MC (the EEX3 option in the KK MC), as 
a function of the cut v msx . It is divided by the Born cross-section. The comparisons are the 
same as in Fig. |21 
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Figure 4: Next-to-leading-log contribution (5\ 
of the /C/C MC), as a function of the cut v 
comparisons are the same as in Fig.|21 
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Figure 5: Sub-NLL contribution /3f J -#n L l for tne YFS3ff MC (the EEX3 option of the /C/C 



MC), as a function of the cut v n 
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Figure 6: Sub-NLL contribution j3{ — /3i NLL with an expanded scale to show clearly the 
differences in the NNLL results for cuts up to 0.95. 
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